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ABSTRACT 


Part  One:  A  MODIFIED  LINEAR  PYROLYSIS  EXPERIMENT 

A  linear  pyrolysis  experiment  utilizing  a  heated  porous  plate  is  proposed  and  analyzed. 
Simple  and  chainlike  surface  gasification  processes  are  defined  and  discussed.  The 
conditions  under  which  the  surface  process  is  a  rate  process  or  one  of  near  equilibrium 
are  determined  in  the  general  case  and  for  chainlike  processes.  It  is  shown  that 
pyrolysis  measurements  yield  the  forward  gasification  rate  in  the  limit  of  a  surface 
rate  process  and  the  equilibrium  vapor  pressure  in  the  limit  of  surface  equilibrium. 
The  pyrolysis  rate  of  potassium  chloride  is  calculated. 

Part  Two:  REMARKS  ON  THE  INTERPRETATION  OF  PYROLYSIS  EXPERIMENTS 
ON  AMMONIA  SALTS 

Previous  interpretations  other  investigators  of  linear  pyrolysis  rate  measurements 
for  NH^Cl,  NH^NO^  and  NH^CIO^  assumed  that  the  surface  process  is  an  unopposed 
rate  process.  This  report  presents  an  alternative  interpretation  of  pyrolysis  experi¬ 
ments  which  may  be  valid  even  when  the  surface  process  is  near  equilibrium.  Some 
features  of  a  previously  given  analysis  of  a  modified  pyrolysis  experiment  (Part  One) 
are  employed  to  demonstrate  that  the  difference  between  the  plate  temperature  and 
the  surface  temperature  of  the  solid  may  be  large  enough  to  affect  seriously  ttie  aotiva- 
tlon  energies  which  were  calculated  from  the  orevious  linear  pyrolysis  data. 
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NOMENCLATURE 


Cj  y  (RgT)/(  2irW.)  ,  average  velocity  of  molecules  of  species  1  In  the  -x 

direction  (cm/sec) 

c  .  specific  heat  at  constant  pressure  of  gaseous  species  i  (cal/g’K) 

N 

Cp  (  ai/aT)Y  N)  ^  ‘^p,  l^i  '  average  specific  heat  of  the  gas 

^  ’  i=l  ’ 

(cal/g»K) 

2 

D  binary  diffusion  coefficient  (cm  /sec) 

2 

H  h  +  (v  /2)  (R  /R  )  ,  total  stagnation  enthalpy  per  unit  mass  (cal/g) 

1  3 

N 

h  ^  h.Yj  ,  total  enthalpy  per  unit  mass  (cal/g) 

i^t  * 

h^  total  enthalpy  per  unit  mass  of  gaseous  species  i  (cal/g) 

total  enthalpy  per  unit  mass  of  the  solid  (cal/g) 

N 

J  y  Y.  h.(T  )  -  h  (T  )  +  (R,/R_)  (vf/2 -v^/2)  (cal/g);  see  Eq.  (1.4) 

^  1|S  I  S  JT  O  13  8  n 

i=l 

J'  seep.  1-3,  above  Eq.  (1.4) 

I  thickness  of  the  gas  film  (cm) 

M  total  number  of  chemical  species  leaving  the  solid  surface 

m  mass  flow  rate  per  unit  area  (g/cm  sec) 

N  total  number  of  gaseous  chemical  species 

Pj  permeability  of  the  porous  plate  to  species  i  (g/cm  sec  atm) 

p  pressure  (atm) 
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partial  pressure  of  species  i  at  the  downstream  side  of  the  porous  plate 
<atm) 

a  constant  in  the  momentum  equation  (atm);  see  Appendix  A 
universal  gas  constant,  1.986  (cal/mole'^K) 

3 

universal  gas  constant.  82. 06  (atm  cm  /mole'K) 

7  2 

universal  gas  constant,  8.317  x  10  (g(cm/8ec)  /mole*K) 
temperature  (®K) 

thickness  of  the  porous  plate  (cm) 
mass -average  velocity  (cm/sec) 
molecular  weight  of  species  i  (g/mole) 

mass  of  molecules  of  type  i  leaving  a  unit  area  of  surface  per  second 
(g/cm\ec) 

spacial  coordinate  (cm) 


ITA 

^  Yj  ,  mass  fraction  of  material  leaving  the  solid  surface 


mass  fraction  of  gaseous  species  i 

accommodation  coefficient  (dimensionless) 
M 

y  0,Y.  /Y 

ZL  i  I.s  s 


«i(pCi)g  /  [dx/(pD)]  ,  ratio  of  a  characteristic  diffusion  time  to  a 
0 

characteristic  recondensation  time 
mass  flux  fraction  of  species  i 

(4/3)p  +  X  ,  effective  viscosity  coefficient  (g/cm  sec) 
coefficient  of  bulk  viscosity  (g/cm  sec) 
thermal  conductivity  of  the  gas  (cal /cm  sec'K) 
ordinary  (shear)  viscosity  coefficient  (g/cm  sec) 
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m  j  [dx/diD)^  ,  dimensionless  spacial  coordinate 
0 


J^[dx/(pD)] 


density  (g/cm  ) 


,  dimensionless  mass  flow  rate 


^  cr^  e.  ,  a  dimensionless  measure  of  the  pressure  at  the  porous  plate 
i=l  ‘ 


I  (PjPp/t)  j  [dx/(pD)]  I  ,  dimensionless  permeability  with  respect 
0 

to  species  i 

dimensionless  surface  temperature;  see  Eq.  (1.15) 

N 

^  0^0. /a  ,  dimensionless  pressure  on  the  downstream  side  of  the  porous 

^1 


(PiPi  J  [<ix/(pD)J  ,  dimensionless  partial  pressure  of  species 


on  the  downstream  side  of  the  plate 


m 

^  ,  dimensionless  forward  surface  gasification  rate 

i=l 

Wi/[a(i(pCi)s  ]  ,  dimensionless  forward  surface  gasification  rate  of 
species  i 


for  surface  equilibrium 
a  gaseous  chemical  species 
a  gaseous  chemical  species 
at  the  upstream  boundary  of  the  solid  (x  » 
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p  at  the  upstream  edge  of  the  porous  hot  plate 

r  a  reference  condition  at  which  the  equilibrium  vapor  pressure  is  known 

8  at  the  solid-gas  interface 

w  for  the  solid 

«o  at  the  downstream  edge  of  the  porous  hot  plate 
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Part  One 

A  MODIFIED  LINEAR  PYROLYSIS  EXPERIMENT 

1.  INTRODUCTION 

Schultz  et.  al.  (Refs.  1-6)  have  reported  many  experiments  on  the  linear  regression 
rates  of  solids.  These  experiments  have  yielded  pyrolysis  rate  data  for  pressures 
and  temperatures  approaching  (but  not  yet  equalling)  those  encountered  in  solid  pro¬ 
pellant  rockets.  These  data  are  of  considerable  practical  importance  in  predicting  the 
deflagration  rate  of  ammonium  perchlorate  (Ref.  7,  pp.  612-618;  Refs.  8-11)  and, 
presumably,  the  burning  rates  of  solid  propellants  (Ref.  12).  Although  the  pressure 
is  sufficiently  high  in  these  experiments  so  that  the  continuum  equations  of  fluid  me¬ 
chanics  may  be  applied  to  the  gas  adjacent  to  the  solid,  the  geometry  of  the  gas  flow 
is  so  complicated  that  unquestionable  Interpretations  of  the  significance  of  the  meas¬ 
urements  appear  to  be  difficult  to  obtain.  Consequently,  it  is  of  Interest  to  investigate 
the  possibility  of  modifying  these  experiments  in  order  to  yield  data  which  can  be 
interpreted  more  easily. 

In  Part  One,  an  analysis  is  presented  for  a  modified  linear  pyrolysis  experiment  which 
should  not  be  difficult  to  perform.  Conditions  under  which  the  surface  process  is  a 
rate  process  or  one  of  near  equilibrium  are  determined  for  the  modified  experiment, 
and  interpretations  of  the  significance  of  the  measurements  are  given  for  these  two 
limiting  cases.  The  pyrolysis  rate  of  potassium  chloride  is  computed  as  an  example. 
The  mathematical  simplicity  of  the  system  described  here  provides  a  strong  incentive 
to  construct  such  an  apparatus  for  measuring  regression  rates.  The  similarity  of  the 
present  model  to  the  previous  pyrolysis  experiments  also  enables  us  to  draw  some 
conclusions  concerning  the  previous  measurements.  * 


•See  Part  Two. 
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2.  MODIFIED  EXPERIMENT 


It  is  suggested  that  the  hot  plate  used  in  previous  pyrolysis  experiments  be  replaced 
by  a  porous  hot  plate.  If  suitable  shielding  can  be  provided  to  prevent  gas  from  es¬ 
caping  without  passing  through  the  porous  plate,  then  the  flow  would  be  approximately 
one-dimensional,  and  the  arrangement  would  appear  as  Illustrated  in  Fig.  1-1. 
Material  gasifies  from  the  solid  surface  at  x  =  0  and  flows  through  the  porous  plate 
at  X  =  /  .  In  analyzing  this  experiment,  the  gas  phase  and  the  interface  conditions 
will  be  considered  separately. 


OS  p  «e 


Fig.  1-1  Schematic  Diagram  of  Modified  Pyrolysis  Experiment 
3.  ANALYSIS  OF  THE  GASEOUS  REGION 


The  one-dimensional,  steady-state  ideal  gas  conservation  equations  for  total  mass, 
the  mass  of  species  i  ,  and  energy  are,  respectively,  * 


pv  =  m  =  constant 


(1.1) 

(1.2) 


•The  basic  equations  are  given  In  Ref.  13,  h).  234-241,  for  exaiiq)le.  Symbols  are 
defined  in  the  Nomenclature.  The  momentum  equation  is  considered  In  Appendix  A. 
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M.  -± 

dx  dx  y  dx  y 


In  Eq.  (1.2),  chemical  reactions  In  the  gas  region  have  been  neglected,  and  the  binary 
diffusion  coefficients  of  all  pairs  of  chemical  species  have  been  assumed  to  be  equal. 
Equation  (1.3)  relies  upon  the  latter  assumption,  the  approximation  that  the  Lewis 
number  Is  unity  (X  =  pD  c^)  ,  and  the  assumption  that  either  the  mass  dlffuslvlty 
equals  the  effective  coefficient  of  viscosity  (tj  =  pD)  or  the  kinetic  energy  is  negli¬ 


gible  in  comparison  to  the  thermal  plus  chemical  energy 


-■'■I- ')  ■ 


A  boundary  condition  for  Eq.  (1.3)  at  x  =  0  is  the  requirement  that  enough  heat  must 

be  conducted  into  the  solid  or  generated  at  the  solid  surface  by  viscous  dissipation  to 

raise  Its  temperature  from  (at  x  = )  to  (at  x  =  0)  and  cause  gasification 

(including  the  accompanying  enthalpy  and  kinetic  energy  changes)  at  T  ,  l.e. , 

8 

(MT/dx)^  -  "I'T,)  -  ■»!>,  (V 


*  ■"  E  2  (’s  ■  (dv/dx)] 


s  mJ' 


which  may  be  written  in  the  form 


(^/Cp)ldH/dx))  ^  -  m  -  "W  *  (“i/Ka 


»  mJ 


♦It  is  to  be  noted  that  c  ,  X  ,  and  p  D  may  depend  upon  both  T  and  Y,  In  this 
formulation.  ^  1 
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where  the  last  term  In  the  curly  brackets  is  usually  negligible.  The  boundary  condi¬ 
tions  for  Eq.  (1.2)  at  x  =  0  are 

=  m£j  ,  i  «  1  ,...,  N  (1*6) 


<s 


where  is  the  mass  flux  fraction  of  species  i  at  the  surface.  *  Employing  Eqs.  (1.4) 

and  (1.6)  and  denoting  the  boundary  value  at  x  =  i  as  Y.  =  Y.  and  H  =  H  , 

1  I*  P  P 

one  finds  the  solutions  to  Eqs.  (1.2)  and  (1.3)  to  be 


(1.6) 


and 

H  =  Hp  -  j(e^P  -  e^) 

X  X 

where  I  =  mf  ^dx/^D)]  and  4  =  mf  [dx/(pD)] 

0  ^0 


(1.7) 


4,  CONDITIONS  AT  THE  POROUS  PLATE 


While  the  experimenter  often  has  fairly  direct  control  over  the  plate  temperature  Tp  , 

the  values  of  Y,  f  which  are  needed  along  with  T  to  determine  H  and  Y. 

i,p  *-  P  P  liP 

In  Eqs.  (1.  6)  and  (1.7)]**  will  basically  be  determined  by  a  permeability  of  the  porous 
plate  to  species  1  ,  P^  ,  and  the  partial  pressure  of  species  1  on  the  downstream 
side  of  the  plate,  p,  .  A  reasonable  phenomenological  expression  is  obtained  by 


♦The  values  are  determined  by  the  overall  stoichiometry  of  the  gasification 
process;  ^  =  0  for  inert  gaseous  species  (i.  e. ,  those  not  leaving  the  solid 
surface). 

•♦Strictly  speaking,  either  Vg  must  also  be  measured  or  the  momentum  equation 
must  also  be  solved  (Appenmx  A)  in  order  to  determine  completely  the  constants 
in  Eq.  (1.7).  However,  the  kinetic  energy  will  usually  be  negligibly  small. 
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assuming  that  the  mass  flux  is  proportional  to  the  pressure  difference  across  the 
plate  for  each  species,  viz. . 


-  (p,/t) 


i.  « 


.1  =  1 . N 


(1.8) 


where  p^  is  the  total  pressure  at  x  =  i,  t  is  the  thickness  of  the  porous  plate, 
and  the  first  term  inside  the  curly  brackets  is  clearly  the  partial  pressure  of  species 
1  at  X  =  1.  In  view  of  Eqs.  (1. 2)  and  (1. 5),  Eq.  (1.  8)  may  be  written  in  the  dimen¬ 
sionless  form 


=  +  CT.i/).  ,  i  =  1 . N  (1.9) 


where  * 


Although  only  N-1  of  the  Y,  are  linearly  independent,  there  are  N  independent 

If  P 

relations  in  Eq.  (1.9)  since  summing  Eq.  (1.9)  over  all  species  i  gives  the  nontrivial 
equation 


where 


«^(l  p  )  =  1 


(1.10) 


N 

O’*  S  o.c.  and  ip 

I  =  1  »  1 


N  N  N 

=  E  o  E  Pi«,/  E 

1  =  1^*  i  =  1  i  =  1 


I 

[Pj[/t^l‘)]  /  Cdx/(pD)] 
0 
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This  result  is  a  consequence  of  the  fact  that  Eq.  (1. 8)  determines  in  terms  of 
N 

p,  and  properties  of  the  porous  plate.  If  Eq.  (1. 10)  is  used  to  compute 

1  =  1 

the  dimensionless  plate  pressure  a  .  then  the  remaining  N-1  Independent  relations 

in  Eq.  (1. 9)  will  determine  Y. 

i»  p 

5.  CONDITIONS  AT  THE  SOLID  SURFACE 

A  general  surface  condition,  obtained  by  equating  the  mass  flow  rate  of  each  species 
to  the  difference  between  the  mass  of  that  species  leaving  and  entering  the  surface 
per  second  (Ref.  14;  Ref.  15,  pp.  60-71;  Refs.  16  and  17)  is 

m€j  =  w.  -  ttjY.  s(pcj)g  ,  i  =  1 . N  (1.11) 

where  w^^  is  the  mass  of  molecules  of  type  i  leaving  a  unit  area  of  surface  per 
second,  is  a  surface  accommodation  coefficient  for  species  i  ,  Cj  =  ^(R^T)/(2irVf^) 
is  the  average  velocity  of  molecules  of  species  i  (of  the  ideal  gas  mixture)  in  the 

-X  direction,  and  p  =  p/  [RgT  ^Zj^(Yj/Wj)]  according  to  the  ideal  gas  law.  Al¬ 
though  the  N  independent  relations  in  Eq.  (1. 11)  should  determine  m  and  the  N-1 
independent  €.  ,  in  most  practical  cases  the  overall  surface  stoichiometry  is  known 
more  accurately  than  w^  or  o:.;  therefore,  Eq.  (1. 11)  would  be  used  to  determine 
m  and  to  provide  N-1  relations  between  N-1  of  the  w^  and  N-1  of  the 
Equation  (1. 11)  may  be  written  in  the  dimensionless  form 

^p^  =  -  Yj  g)  .  i  =  1 . N  (1.12) 


♦It  should  be  emphasized  that  wj  and  oj  may  be  functions  of  Tg  ,  Pg  and  all  Y|  g 
for  complex  surface  processes.  ’ 
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I 

where  5  aj{pCj)g  /  [dx/(^)]  and  Wj  =  w^/ ]  .  From  Eqs.  (1.11) 

0 

or  (1. 12)  It  follows  that.  In  general,  the  surface  process  with  respect  to  species  i 
will  be  a  rate  process  similar  to  that  encountered  In  vacuum  vaporization  or  sublima¬ 
tion  experiments  (Refs.  17  — 19)  if 

while  near  equilibrium  conditions  prevail  for  species  1  at  the  surface  If  ~  1» 

that  is ,  If 

(4p€j)/(Wi^j)  «1  . 

6.  GENERAL  EQUATIONS  FOR  DETERMINATION  OF  PYROLYSIS  RATE 

If  Eqs,  (1.6)  and  (1.7)  are  evaluated  at  x  =  0  ,  and  if  Eq.  (1.6)  is  solved  for  Y. 

l»  ® 

and  Eq.  (1.7)  solved  for  i  ^  ,  then  it  is  found  that 

\b  =  '^l.p®  ®  ^)  •  1  =  1 . N  (1.13) 

and 

T  =  e  -  1  (1.14) 

where 

T  “  (Hp  -  Hg)/J  (1. 15) 

Utilizing  Eq.  (1. 14),  one  may  also  write  Eq.  (1.  13)  in  the  form 

^I,p  “  ^1,8  ""  ^^1.8  ~  ^j) i  =  1  . . .  • .  N  (1.13a) 


1-7 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


3-64-81-1 


If  Eq.  (1. 13a)  Is  multiplied  by  h,(T  )  and  summed  over  all  i  and  if  use  is  made  of 

1  s 

Eq.  (1. 15)  and  the  definitions  of  H  and  J,  then  by  solving  the  resulting  equation  for 

=  T  •  it  is  found  that 


N 

-  fj 

-  (1.15a) 


From  Eq.  (1. 15a)  it  is  seen  that  t  is  physically  a  dimensionless  measure  of  the 

unknown  surface  temperature  T  .  Consider,  for  example,  the  special  case  in  which 

s 

kinetic  energy  and  viscous  dissipation  are  negligible  and  the  specific  heats  of  all  gaseous 

species  are  constant  and  equal,  whence  Eq.  (1. 15a)  becomes  t  =  (T  -  T  )/(J'/c  )  . 

p  s  p 

Equations  (1.9),  (1.12),  (1.13)  and  (1.14),  plus  the  tvyo  equations  obtained  by 
summing  Y  over  i  at  the  surface  and  at  the  plate,  basically  constitute  3N  +  3  independ¬ 
ent  equations  in  the  3N  +  3  unknowns  Y,  ,  a  ,  i  ,  e.  ,  Y.  ,  and  t  .  Appear- 

1  >  p  pi  i » s 

ing  in  these  equations  are  the  quantities  ,  ip^  ,  tp  ,  0^  ,  and  ,  which 

involve  known  parameters  and  known  functions  of  the  unknown  variables.  Hence  it  is 
presumably  possible  to  solve  this  set  of  equations  for  the  unknown  variables  and,  in 
particular,  to  compute  the  dimensionless  pyrolysis  rate  ^  p-  *  Essential  features  of 
the  solution  will  best  be  illustrated  by  considering  special  simple  cases. 


*In  addition  to  the  unknown  pyrolysis  rate,  $  actually  also  contains  the  length  £ 
which  may  be  unknown  in  some  experiments.*^  In  such  cases  the  momentum 
equation  eventually  determines  £  (Appendix  A). 
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7.  SIMPLE  CHAINLIKE  SURFACE  PROCESSES 

Many  of  the  Important  characteristics  of  this  problem  may  be  deduced  by  examining 
the  simplified  set  of  equations  henceforth  referred  to  us  the  set  (I)  , 


»♦  + V 


1  -  (I  -  Vp)/(1  +T) 


=  ln(l  +  T) 


Ip  =  -  Y^) 


(1.16) 


(1.17) 


(1.18) 


(1. 19) 


(1.20) 


where  j/i  is  a  known  constant,  Y  =  Y,  ,  is  the  total  mass  fraction  of  all 

i  =  1  ^ 

chemical  species  with  Cj  (while  =  o  for  i  =  M  +  1,  .  .  .  ,  N)  ,  and 

M  M  /  M  \  i 

(3=  j  ^i^i,  s'^^s  w  H  Z)  [dx/(pD)  3  /j9  are  known 

functions  of  t  and  a  alone.  The  general  equations  reduce  to  the  set  (I)  for  the 
simple  chainlike  surface  processes  discussed  below. 

If  parallel  competing  surface  reactions  do  not  occur,  a  condition  under  which  the  sur¬ 
face  process  will  be  called  chainlike,  then  the  single  concentration  variable  Y  (which 
includes  all  of  the  chemical  species  leaving  the  solid  surface)  may  be  sufficient  to 
characterize  the  composition  of  the  gas  for  purposes  of  pyrolysis  rate  calculations. 

M 

Since  e  =  1  ,  Eq.  (1.18)  can  be  obtained  by  summing  Eq.  (1.13)  from  1  =  1 
1  =  1^ 
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to  i  =  M  and  utilizing  Eq.  (1.14)  to  eliminate  |p  .  Formally,  Eq.  (1.20)  is  obtained 
by  summing  Eq.  (1. 12)  from  i  =  1  to  1  =  M  .  Equation  (1, 17)  can  be  derived  by 

summing  Eq.  (1.9)  from  i  =  1  to  i  =  M  when  =  Wj^(i  =  1 . N)  , 

P ,  =  P-(i  =  1 . N)  ,  and  p  =  0(i  =  1 . M);  if  any  of  these  conditions 

are  violated,  an  expression  for  which  is  somewhat  more  complicated  than  Eq.  (1. 17) 
is  obtained.  The  two  remaining  relations  in  set  (I)  ,  Eqs.  (1. 16)  and  (1. 19)  are  merely 
alternative  forms  of  Eqs.  (1.10)  and  (1.14),  respectively. 

Further  restrictions  are  necessary  in  order  to  satisfy  the  additional  requirement  that 
p  and  oj  must  depend  only  upon  t  and  a  .  The  surface  process  may  be  called 
simple  if  w,  and  a.  depend  upon  T_  only,  since  gaseous  concentrations  (p_  ,  Y,  _) 

1 _ 1 _ S _  S  1|  B 

do  not  affect  the  surface  rates  in  this  case.  Assuming  =  W^(i  =  1  , . . .  ,  N)  and 
i 

holding  J  Ldx/(pD)J  constant,  it  follows  from  the  definition  of  and  oj.  that 
0  ^ 

"i  ^  "i^'^s^  implies  =  ^.(T^  ,  p^)  and  that  =  Q;^(Tg)  and  w^  =  Wj(Tg) 

together  imply  ,  p^)  ,  whence  the  additional  assumption  =  p(i  =  1 . M) 

implies  that  /3  =  /3  (T  ,  p  )  and  w  =  w(T  ,  p  )  according  to  the  definitions  of  /3 
s  s  s  s 

and  o)  .  If  one  neglects  kinetic  energy  and  viscous  dissipation  and  assumes  that  the 
heat  capacities  of  all  species  (i  =  1  , . . . ,  N)  equal  c^  between  T^  and  Tp  ,  it 
follows  from  Eq.  (1. 15a)  that  (with  T  given)  T  =  T  (t)  .  The  momentum  equation 

S  5 

(Appendix  A)  implies  that  p  is  determined  by  p  ,  and,  therefore,  from  the  definition 

s  p 

of  (7  .  P_  =  p„(a)  .  Hence,  under  the  conditions  stated  above,  B  =p  {t,  a)  and 
s  s 

W  =  w(t,  o)  .  * 


•Some  consequences  of  the  approximations  Pi  =/S(i  =  1,...,  M)  ,Wi  =  Wi(i  =  1, . . . ,  N), 

Pi  =  Pi(i  =  1 . N)  ,  and  pi^«,  =  0(i  =  1 . M)  are  Yi  p  =  ^  Yp(i  =  1 . M) 

Csee  Eqs.  (1,9)  and  (1.17)],  Yi^g  =  ^  Ys(i  =  1 . M)  Esee’Eqs.  (1.13)  and  (1.18)3, 

and  fa)i  =  =  1, . . M)  [see  Eqs.  (1.12)  and  (1.20)3  .  While  the  present  approxi¬ 

mations  are  not  unreasonable,  in  some  cases  it  is  possible  to  obtain  u  =  cu(T,a)  and 
P  =  p{T,  df  as  well  as  the  set  (I)  under  much  less  stringent  conditions. 
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8.  THE  FUNCTIONS /3  (<r.  t)  and  oKff.T)  FOR  SIMPLE  CHAINLIKE  SURFACE 
PROCESSES 

For  simple  chainlike  surface  processes  it  is  actually  possible  to  obtain  the  T  and 

B 

p  (i.e.,  T  and  a  )  dependence  of  a>  explicitly.  In  general,  in  Eq.  (1.11) 
m  =  m(T  ,  Y.  „  ,  p„)  ,  whence,  at  constant  p„  ,  Eq.  (1.12)  becomes 

8  8  S  S 

‘  1 . « 

for  simple  surface  processes.  But,  in  surface  equilibrium,  which  is  defined  by 

^p  ^  0  ,  it  is  a  consequence  of  thermodynamics  that,  for  chalnllke  surface  processes, 

Y.  =  Y.  ^(T„)  ,  where  the  subscript  e  identifies  equilibrium  conditions  at  the 

1)8  1|S|6  S 

given  constant  p  .  Hence, 
s 


€.  Urn  [|p(Tg  .  Yj  g)]  =  0  i,  j  =  1 . M  (1.22) 


Y  —  Y 
J.s  j,s,e 


T  fixed 
s 


Substituting  Eq.  (1.21)  into  Eq.  (1.22),  one  gets 


‘^i^V  =  Yi,s,e(V  •  '  =  1 . ^ 


(1.23) 


Summing  Eq.  (1.23)  from  i  =  1  to  i  =  M,  one  gets 


M 


w(T 


s)  =  Yi,3,e(Tg)  =  Yg  g(Tg) 


(1.24) 
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when  /3j  a  j9  .  Assuming  that  the  gas  Is  Ideal,  the  density  of  the  gas  Is  negligible 
compared  to  that  of  the  solid,  and  Wj  =  Wj^(l  =  1  , . . . ,  N)  ,  the  Clapeyron  equation, 
becomes  (Ref.  20 ,  pp.  450-460) 


Hs,A,e»]  =  *l[,E  -  \(v]dV(RjT^)  (l.i 


for  this  surface  process.  Integrating  Eq.  (1.25)  at  constant  P„  „  =  P  (  =  P_  J  from 
a  temperature  T  ,  at  which  the  vapor  pressure  (Y  p  .)  is  known,  to  temperature 

I*  0  0 

T  ,  one  gets 
s 


Equations  (1. 24)  and  (1.26)  give  the  T  dependence  of  o)  explicitly.  There  is  some 

s 

evidence  that  p  is  not  strongly  dependent  upon  T  at  constant  p  for  a  few  simple 

s  s 

chainlike  surface  processes  (Ref.  15,  pp.  320-324).  It  is  clear  from  the  definition 

of  /5j  and  that  ~Pg  constant 

T  ,  whence  also  /3  ~  p  and  w  ~  1/P„  • 
s  s  s 

9.  IMPLICATIONS  OF  THE  SET  (I)  REGARDING  THE  GASIFICATION  PROCESS 

Using  Eq.  (1.16)  for  ct  in  Eq.  (1.17),  the  result  for  Yp  in  Eq.  (1.18),  this  result 

for  Y  in  Eq.  (1.20),  and  Eq.  (1. 19)  to  eliminate  |  from  the  final  expression, 
s  p 

it  is  found  that 

ln(l  +T  )  =/3  I  w  -  1  +  (1  +  T)"^  [l  +  {l/>f>)  ln(l  +  T)]"^|  (1.27) 

which  may  be  solved  for  t  ,  after  which  Eqs.  (1. 16)  through  (1. 19)  give  the  remain¬ 
ing  variables. 
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In  the  limiting  case  t«  1  •  *  for  all  values  of  ip  except  ^  «1  ,  Eq.  (1.19)  implies 


Ip-  r  («1) 


(1.28) 


Eqs.  (1. 16)  and  (1. 17)  imply 


Yp  -  t/ 4!  {«!) 


(1.29) 


and  Eq.  (1.27)  becomes** 


T  «  00)7  1  +  /3  [l  +  (1/i/O] 


(1.30) 


Since  the  conditions  for  near  equilibrium  surface  processes  and  unopposed  surface 
rate  processes  given  in  Section  5  reduce  to 


|p/(co/3)  «  1 


1  -  [l  ]  «  1 

respectively,  for  chainlike  surface  processes,  it  is  seen  from  Eqs.  (1.28)  and  (1.30) 
that  when  t  «  1  the  quantity 


|p/(aj/3)  =  I  1  +  /?[l  +  (1/i/;)] 


(1.31) 


*Thls  has  been  true  in  many  of  the  pyrolysis  experiments  reported  previously. 
** From  Eqs.  (1.18),  (1. 20)  and  (1. 30),  it  is  easily  shown  that  PI  +  (l/fl)  + 
(1/*^)]  »  cu  >  Y«  >  Yp  whL  0  <  T  «  1  . 
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determines  whether  equilibrium  or  rate  conditions  are  valid  at  the  surface.  From 
Eq.  (1..31)  it  is  seen  that  (since  '/'  is  not  1)  near  equilibrium  conditions  will 
prevail  at  the  surface  only  i!"  p  »  \  ,  while  the  surface  process  will  be  an  unop¬ 
posed  rate  process  only  if  P  «  1  .  From  the  definition  of  p  ,  which  represents 
the  ratio  of  a  characteristic  diffusion  time  to  a  characteristic  recondensation  time, 
it  can  be  seen  that  high  pressure,  large  a  ,  a  large  accommodation  coefficient,  and 
low  diffusivities  favor  surface  equilibrium.  * 

In  most  pyrolysis  experiments  (Refs.  1  —  6),  p  =  1  atm  and  300“ K  %  900“ K. 

4  s 

For  reasonable  pD  and  (pD  ^  5  x  lO  g/cm  sec,  W  «  75  g/ mole), 

/  4  -1  ^  ^ 

P/a  2  =  5  XIO  cm  within  a  factor  of  10.  Since  it  is  expected  that  one  mean  free 

-5  -5  -2 

path  ~  10  cm,  it  is  therefore  estimated  that  10  cm  ^  i  ^10  cm,  and  it 

1  3 

follows  that  10  0!  ~  ^  ~10  or.  Hence,  for  simple  chainlike  surface  processes 

with  T  «  1  ,  the  surface  process  will  probably  be  near  equilibrium  for  a  =  1  and 
a  rate  process  for  o;  <  0.01  . 

10.  INTERPRETATION  OF  PYROLYSIS  RATE  MEASUREMENTS 

When  the  surface  rate  controls  the  pyrolysis  rate  (Section  9),  then 

4p  =  oop  (1.32) 

M 

which  implies  m  =  ^  w,  according  to  the  definition  of  w  and  |  .  Hence  pyrolysis 

i  =  1  ^  P 

rate  measurements  yield  directly  the  forward  surface  gasification  kinetics.  This  has 

been  the  basis  of  interpretations  of  previous  pyrolysis  data  (Refs.  1,3,  5,  21  —  23). 

♦Since  use  has  not  been  made  here  of  the  assumption  that  the  surface  process  is 
simple,  these  results  are  valid  for  all  chainlike  gasification  processes  (with 
T  «  1  and  ^  not  «  1). 
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When  the  surface  process  Is  near  equilibrium,  the  left-hand  side  of  Eq.  (1.201  is 
much  less  than  either  term  on  the  right;  therefore  the  left-hand  side  of  Eq.  (1.27)  is 
negligible.  Utilizing  Eq.  (1.19)  in  Eq.  (1.27),  one  gets 


+  U 


1  -  a; 


(1.33) 


which  shows  that  (knowing  ip  )  a  measurement  4p  determines  u)  .  When  t  «  1 
(and  jp  not  «1),  Eq.  (1. 28)  Implies  that  Eq.  (1. 33)  may  be  solved  explicitly  for 
by  expanding  the  left-hand  side  in  powers  of  ^  p  8-od  neglecting  terms  of  higher  than 
first  order,  viz. , 


Ip  «  oo/ll  +  (1/0)] 


(1.33  a) 


Equation  (1.33a)  is  an  expression  which  may  also  be  derived  from  Eqs.  (1. 28)  and 
(1.30),  without  using  Eq.  (1.33),  by  utilizing  the  previous  result  that  ^  »  1  near 
surface  equilibrium  when  t  «  1.  Since  oj  =  Y  at  surface  equilibrium  fEq.  (1.20)], 
it  follows  from  Eqs.  (1.33)  and  (1.26)  that  for  near  equilibrium  surface  processes  a 
pyrolysis  rate  measurement  yields  the  equilibrium  vapor  pressure  (provided  i  ,  pD  , 
and  ip  are  known),  and  that  the  dependence  of  the  pyrolysis  rate  upon  determines 
the  heat  of  gasification  (even  if  i  ,  p  D  ,  and  tp  are  not  known,  so  long  as  they  remain 
constant) . 


11.  PYROLYSIS  RATE  OF  POTASSIUM  CHLORIDE 


As  an  example  of  the  application  of  the  preceding  analysis  we  shall  consider  KCl  , 
one  of  the  few  materials  whose  properties  are  sufficiently  well  known  to  permit  an 
a  priori  calculation  of  the  pyrolysis  rate.  Sublimation  of  KCl  is  a  slnq)le  chainlike 
surface  process  (Ref.  22).  In  Fig.  1-2,  the  theoretical  pyrolysis  rate  is  plotted  as 
a  function  of  the  plate  temperature  T^  for  various  values  of  i  as  con^nited  from 
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the  set  (I)  [using  specifically  Eqs.  (1.19),  (1.24).  (1.26).  and  (1. 27)]wlth  1/0  =  0. 

T  =  300*K.  and  p  =1  atm  .  Here  (Refs.  22.  24.  and  25)  c  =  0.117  cal/gTC, 
o  *^s  P 

N 

h^(Tg)  -  h^(T^)  =  0.166(Tg  -  T^)  cal/g  .  '  W  "  * 

pD  =  5  X  10  ^  g/cm  sec.  a  =  0.63  ,  and  Y  =  6.4  x  10^  exp(-26.070/T  )  . 
The  computations  were  completed  only  for  T  less  than  the  melting  point  (1049*  K). 

S  _g 

The  results  imply  that  surface  equilibrium  exists  (/3  >  1  provided  £,  >10  cm) 
and  show  that  (Y  -  Y  )/Y  <  0.7  (Fig.  1-3).  Since  it  was  always  found  that 

S I  @  S  S  •  0 

T  -  T  <  5*K,  Eq.  (1. 33a)  is  actually  valid  here,  and  a  pyrolysis  rate  measure- 
P  s 

ment  yields  the  equilibrium  sublimation  pressure. 
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Appendix  A 

THE  MOMENTUM  EQUATION 

In  the  gaseous  region,  the  one-dimensional  momentum  equation  is 


_  A. 


Rgy  dx  dx 


(-i)= « 


(A.l) 


the  first  Integral  of  which  is 


(R-/R„)p  +  mv  ~  mdv/d4  =  (R  /R  )p  =  constant 


(A.  2) 


when  Tj  =  pD  .  The  evaluation  of  the  constant  p  depends  upon  the  experimental 

IT 

arrangement.  If  the  solid  is  supported  in  such  a  way  that  the  spacing  i  is  fixed,  * 
then  p  must  be  evaluated  in  the  gas  (e.  g. ,  at  x  =  i  )  .  If,  however,  only  the  gas 

TT 

layer  supports  the  solid  (i.e. ,  i  is  maintained  by  the  pressure  in  the  gas  layer),  then 
the  momentum  equation  extends  into  the  solid  and  p  may  be  evaluated  in  the  solid, 

TT 

where  velocity  gradients  are  negligibly  small  (i.e. ,  viscous  effects  are  negligible) 

3 

and  the  dynamical  pressure  (mv  )  is  negligible  (since  p  ^  10  p  implies 
g  ^  TT  gas 

V  ~  10  Thus,  in  the  latter  case,  p  is  the  hydrostatic  pressure  in  the 

ir  go-s  7r 

solid,  which  equals  the  atmospheric  pressure  plus  any  additional  pressure  (applied 
force  divided  by  cross-sectional  area  of  sample)  applied  to  the  solid  sample  in  order 
to  increase  "thermal  contact"  with  the  porous  hot  plate.  While  Eq.  (A.  2)  and  the  ideal 
gas  equation  of  state  comprise  the  additional  relations  required  to  determine  the  dis¬ 
tribution  of  all  other  properties  (in  particular,  p  and  v  )  in  the  gas,  Eq.  (A.  2) 
cannot  be  integrated  analytically. 


♦This  appears  to  be  the  case  in  the  previous  pyrolysis  experiments. 
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Considering  the  case  In  which  viscous  dissipation  in  the  gas  is  negligible,  Eq.  (A.  2) 
explicitly  gives  v  In  terms  of  p  ,  viz. , 


V 


^2  ^ 


-  p)/m 


Utilizing  Eq.  (1. 1)  and  the  ideal  gas  law  in  Eq.  (A-3),  we  find 


(A.  3) 


(A.  4) 


which  determines  p  in  terms  of  the  variables  Y^.  and  T  .  Equations  (A.  3)  and  (A.  4) 

may,  of  course,  be  used  to  evaluate  >  Pp  >  Pg  terms  of  p  ,  , 

Y.  ,  ,  m  ,  and  p  . 

i,  s  s  tt 

It  may  be  noted  that  Eq.  (1. 10)  and  the  relation  obtained  by  evaluating  Eq.  (A.  4)  at 

X  -  I  constitute  two  independent  equations  for  the  single  unknown  p^  (i.e. ,  a)  . 

By  equating  the  values  of  Pp  given  by  these  two  equations,  an  expression  is  obtained 

which  may  be  solved  for  p^  if  m  is  known  or  for  m  if  p^  is  known.  In  the  former 

case  m  would  be  determined  by  the  remaining  governing  equations  (through  |  p  ) 

with  I  given,  while  in  the  latter  case  the  pyrolysis  rate  is  determined  by  the  dynamic 

aspects  of  the  problem  and  the  remaining  governing  equations  determine  l  (through  |  ) 

It  2  P 

(Section  4) .  This  is  most  easily  illustrated  in  the  case  mv  «  ^  p  (which  is  usually 
valid),  whence  Eq.  (A. 4)  reduces  to 


P  =  P„  =  P„  =  P  (A.  5) 

»  P  TT 

If  Pj  =  Pj  (i  =  1  , . . . ,  N)  and  (i  =  1  , . . . ,  N)  ,  then  summing  Eq.  (1. 8) 

over  all  i  yields 


m  =  (Pj/t)  (p^  -  p^) 


(A.  6) 
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where  use  has  been  made  of  Eq.  (A.  6)  and 


If  ^  is  fixed  so 


that  the  other  governing  equations  determine  m  ,  then  Eq.  (A.  6)  determines  p  ; 

TT 

if  p^  is  known,  then  Eq.  (A.  6)  determines  ra  directly. 
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Part  Two 

REMARKS  ON  THE  INTERPRETATION  OF  PYROLYSIS 
EXPERIMENTS  ON  AMMONIA  SALTS 


1.  INTRODUCTION 

Schultz  and  co-workers  have  reported  linear  pyrolysis  rate  measurements  for  NH^Cl 
(Refs.  1  and  2),  NH^NO^  (Ref.  3)  ,  and  NH^CIO^  (Ref.  4),  and  have  given  consistent 
interpretations  of  their  results  based  on  absolute  reaction  rate  theory  (Refs.  1,3,  and 
and  5).  These  interpretations  rely  upon  the  assumption  that  the  surface  process  is  an 
unopposed  rate  process  (i.e. ,  surface  conditions  are  far  from  equilibrium).  Since  the 
pressure  in  the  gas  film  between  the  pyrolyzing  solid  sample  and  the  hot  plate  is  prob¬ 
ably  sufficiently  high  for  the  continuum  equations  of  fluid  mechanics  to  be  applicable 
in  most  of  these  experiments,  it  is  conceivable  that  equilibrium  surface  conditions 
may  be  approached  and  that  heat  and  mass  transfer  rates  may  be  important  in  deter¬ 
mining  the  measured  pyrolysis  rates.  One  purpose  of  these  comments  is  to  present 
an  interpretation  of  the  linear  pyrolysis  experiments  which  may  be  valid  even  in  the 
limit  in  which  fluid  processes  control  the  pyrolysis  rate. 

A  detailed  analysis  (including  heat  and  mass  transfer)  of  a  model  of  a  modified  pyrol¬ 
ysis  experiment,  in  which  the  gases  pass  through  a  porous  hot  plate,  has  been  given 
(Part  One).  Since  this  model  has  many  features  in  common  with  the  previously 
reported  pyrolysis  experiments,  certain  of  the  results  in  Part  One  may  be  used  in 
interpreting  pyrolysis  rate  measurements.  In  particular,  the  difference  between  the 
plate  temperature  and  the  surface  temperature  of  the  solid  will  be  discussed.  The 
symbols  used  are  defined  in  the  Nomenclature. 
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2.  ALTERNATIVE  INTERPRETATION  OF  PYROLYSIS  EXPERIMENTS 


It  has  previously  been  assumed  that  in  pyrolysis  <d  ammonia  salt.s  the  overall  surface 
process  is  one  of  dissociation  (Refs.  2,  3  and  4),  viz. , 

NH^Cl(s)  ^  NH3(g)  1-  HCl(g)  AH  - 

40.  3  kcal/mole 

(2.1) 

NH^N03(i)  -  NH3(g)  +  HN03(g)  AH  = 

39.  8  kcal/mole 

(2.2) 

and 

NH^ClO^(s)  ^NH3(g)  +  HClO^(g)  AH  = 

56  kcal/ mole 

(2.3) 

This  is  an  hypothesis  which  is  supported  by  some  evidence  for  dissociation  in  vacuum 
sublimation  experiments  (Refs.  1,  3,  and  5-9).  The  AH  values  listed  here  are 
obtained  mainly  from  equilibrium  vapor  pressure  measurements  (Refs.  6,  7,  8  and  10), 
except  in  the  case  of  NH^CIO^  for  which  AH  must  be  calculated  from  the  heats  of 
formation  since  no  vapor  pressure  measurements  exist. 

In  the  pressure  and  temperature  range  of  the  pyrolysis  experiments,  it  may  be  true 
that  either  (1)  the  equilibrium  composition  is  one  with  a  negligible  amount  of  dissocia¬ 
tion,  or  (2)  dissociation  occurs  in  a  gaseous  process,  the  rate  of  which  is  small  com¬ 
pared  with  the  rate  of  escape  of  the  gas  from  the  film  between  the  plate  and  the  sample. 
In  either  of  these  cases,  the  overall  surface  process  would  be  one  of  ordinary  sublima¬ 
tion,  not  accompanied  by  dissociation.  Since  the  surface  accommodation  coefficients 
a  for  these  processes  are  not  known,  *  it  would  still  be  impossible  to  predict,  from 
an  analysis  similar  to  that  given  in  Part  One,  whether  the  surface  process  is  near 
equilibrium  or  a  rate  process.  However,  if  the  surface  processes  are  simple  chainlike 


♦Since  the  appropriate  accommodation  coefficients  here  are  those  for  the  undissociated 
vapor,  previously  reported  accommodation  coefficients  (e.g. ,  Ref.  6)  are  not  applicable. 
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processes  such  as  that  encountered  in  KCl  sublimation  {Part  One) ,  then  it  may  be 
expected  that  the  temperature  dependence  of  the  pyrolysis  rate  would  determine  the 
heat  of  sublimation  regardless  of  the  importance  of  the  fluid  transfer  processes 
(Part  One).  The  pyrolysis  results  (Refs.  2,  3  and  4)*  are 

m/p  =■  120  exp  r-13.  50Q/(R,T  )]  cm/sec  (629*K  <  T  <  807“K)  (2.4) 

for  NH^Cl  , 

m/p  =  120  exp  [-7100/(R  T  )]  cm/sec  (453*  K  <  T  <  573°  K)  (2.5) 

TT  1  P  P 

for  NH^NOg  ,  and 

m/p^  =  4600  exp  [-20,  000/(R^Tp)]  cm/sec  (753*  K  <  <  903“K)  (2.6) 

for  NH^CIO^  .  These  expressions  would  therefore  imply  that** 


NH^Cl(s)  -NH^CKg) 

AH  =  13.  5  kcal/mole 

(2.7) 

NH4N03(i)  -NH4N03(g) 

AH  =  7.1  kcal/mole 

(2.8) 

NH4C104(s)  NH4C104(g) 

AH  =  20  kcal/mole 

(2.9) 

respectively,  in  the  temperature  ranges  indicated  above.  The  AH  values  given  in 
Eqs.  (2.7)  through  (2.9)  would  actually  be  valid  for  arbitrary  types  of  surface  processes 
if  surface  equilibrium  prevails  in  the  pyrolysis  experiments,  but  must  be  restricted  to 
simple  chainlike  processes  if  the  surface  conditions  are  far  from  equilibrium. 


♦Undefined  notation  is  the  same  as  that  used  in  Part  One. 

♦♦These  AH  values  rest  on  the  assumption  that  Tp  differs  little  from  Ts  ;  in  this 
connection  see  Section  3,  Part  Two,  where  improved  aH  values  are  estimated. 
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We  now  briefly  discuss  each  material  separately,  pointing  out  arguments  for  and 
against  the  processes  given  in  Eqs.  (2.7)  through  (2.9). 

A.  NH.Cl 

4 

Since  the  pyrolysis  data  (Ref.  2)  yield  the  same  activation  energy  as  the  vacuum  subli¬ 
mation  data  (Ref.  6)  and  the  low  pressure  equilibrium  appears  at  present  almost 
certainly  to  involve  entirely  dissociated  NH^Cl  in  the  gas  (Ref,  7) ,  it  follows  that 
the  new  interpretation  will  be  valid  only  if  hypothesis  (2)  of  the  previous  paragraph 
holds  true .  The  large  effect  of  HgO  in  catalyzing  the  decomposition  (Ref.  7)  and 
the  absence  of  H2O  from  the  dried  NH^Cl  samples  used  in  the  pyrolysis  experiments 
(Ref.  1)  make  it  probable  that  the  gas-phase  dissociation  is  slow.  The  agreement 
between  the  data  of  Refs.  2  and  6  indicates  that  the  surface  sublimation  is  probably  a 
rate  process  and  a  <0.01  (Part  One).  The  break  in  the  experimental  pyrolysis 
curve  might  result  from  an  increase  in  the  gas-phase  dissociation  rate.  • 

B.  NH.NO„ 

4  3 

In  pyrolysis,  NH^NOg  melts  before  gasification,  and  a  change  in  the  mechanism  of 
gasification  at  the  melting  point  causes  the  break  in  the  experimental  pyrolysis  curve 
(Ref.  3).  Since  existing  evidence  indicates  that  in  equilibrium  the  gas  is  completely 
dissociated  (Ref.  8),  hypothesis  (2)  seems  more  plausible  than  (1)  for  the  temperature 
range  of  the  pyrolysis  experiments. 


^Although  it  has  been  suggested  that  the  break  in  the  pyrolysis  rate  curve  is  caused 
by  diffusion  effects  becoming  rate  controlling  when  ^e  equilibrium  dissociation 
pressure  becomes  less  than  atmospheric  pressure  (Ref.  2),  this  explanation  is 
questionable  on  the  grounds  that  the  equilibrium  dissociation  pressure  is  always 
below  atmospheric  pressure  for  NH^N03  which  exhibits  no  sharp  decrease  in 
ln(m/p^)  with  increasing  l/T  • 
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C.  NH.CIO. 

4  4 

In  spite  of  the  large  amount  of  recent  work  on  NH^CIO^  ,  *  there  appears  to  be  no 
conclusive  experimental  information  to  eliminate  liypotheses  (1)  or  (2).  Evidence 

k- 

that  some  NH.CIO.  condenses  out  of  the  gas,  which  was  interpreted  as  favoring  the 
4  4 

dissociation  process  given  in  Eq.  (2.3)  (Ref.  9),  is  also  consistent  with  Eq.  (2.9). 

It  has  been  pointed  out  (Ref.  11)  that  since  the  very  strong  acid  HCIO^  would  have  a 
much  greater  tendency  to  form  hydrogen  bonds  than  would  HNO^  or  HCl,  it  should 
favor  the  associated  vapor  NH^CIO^  .  Hence,  the  equilibrium  associated  vapor 
hypothesis  (1)  appears  to  have  considerable  theoretical  justification.  If  the  reaction 
in  Eq.  (2.9)  is  valid,  then  the  pyrolysis  experiments  apparently  provide  the  first 
determination  of  any  thermodynamic  property  of  undissociated  NH^ClO^(g)  . 

3.  DIFFERENCE  BETWEEN  PLATE  TEMPERATURE  AND  SURFACE  TEMPERATURE 

In  the  previous  interpretations  of  the  pyrolysis  experiments  it  has  been  assumed  that 
the  plate  and  surface  temperatures  are  approximately  equal  in  order  to  compute  activa¬ 
tion  energies  for  the  surface  process.  It  appears  reasonable  that  in  the  pyrolysis  ex¬ 
periments  heat  is  transferred  from  the  hot  plate  to  the  solid  surface  across  the  gas 
film  by  a  conduction  process  which  is  essentially  one-dimensional.  If  this  is  the  case, 

then  the  difference  between  the  plate  temperature  T  and  the  surface  temperature  T 

P  8 

can  easily  be  estimated  in  terms  of  the  observed  pyrolysis  rate.  The  appropriate  one¬ 
dimensional  heat  conduction  equation  has,  in  fact,  been  solved  in  Part  One  [^Eqs.  (1. 14) 
and  (1. 15a)]  ,  and  the  result  may  be  written  in  the  form 

/r  N  1^1 

M  ^i¥V  ■  W  =  m/[dx/(pD)]  -1  (2.10) 

♦See,  for  example.  Eighth  International  Symposium  on  Combustion.  Baltimore, 
Williams  and  Wilkins  Co.  (in  press). 
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where  use  has  been  made  of  the  fact  that  the  pyrolysis  rates  are  so  small  that  kinetic 
energy  and  viscous  dissipation  are  negligible.  It  should  be  emphasized  that  Eq.  (2. 10) 
relies  upon  none  of  the  aspects  of  the  model  in  Part  One  except  the  heat  transfer 
equation. 

Since  is  1.527  g/cm^  for  NH^Cl  ,  1.725  g/cm'^  for  NH^NO^  ,  and  1.95  g/cm^ 
for  NH.CIO.  ,  it  follows  from  Eqs.  (2.4)  through  (2.6)  that  m  never  exceeds 
0.35g/cm  sec.  Utilizing  the  reasonable  values  p  D  =  5  x  10  g/cm  sec  and 
,  -4 

■*  -  10  cm,  it  is  seen  that  m  J  [dx/(pD)]  <<  1  ,  whence  the  right-hand  side  of 

0 

Eq.  (2. 10)  may  be  expanded  in  powers  of  this  quantity  and  powers  higher  than  the 
first  will  be  negligible.  With  the  reasonable  estimate  c  =0.2  cal/g  ‘K  ,  Eq.  (2. 10) 
then  becomes 

■  N 

Tp  -  \  ,2  'iW  -  W  -K  (2.11) 

- 

The  quantity  in  the  square  brackets  in  Eq.  (2. 11)  may  be  computed  from  the  relations 

<2-12) 

and 

~  -  V  ^  ^t  %.7r.2<\  -  ^t)  (2.13) 

where  W  is  the  molecular  weight  of  the  salt,  is  a  temperature  at  which  a  phase 
transition  of  the  sample  occurs,  Ah^  is  the  heat  of  transition  per  unit  mass,  and 

^p,  IT,  1  °p,7r,  2  ®^P®cities  of  the  sample  below  and  above  the  transi¬ 

tion  temperature,  respectively.  For  NH^Cl  (Refs.  6,  7,  10  and  12),  W  =  53.5g/mole, 
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=  457. 6"  K  (a  solid  phase  ti-msition),  -  1000  cal/mole,  ^  =  20. 1  cal/ 

mole  *K,  and  c  „W  ^  22.3  cal/mole  "K;  for  NH.NO.,  (Refs.  S,  10,  and  12), 

Pi7r,2  4  .1 

W  =  80  g/mole.  442.  H“  K  (mellin}^),  Ah,W  l.SOO  cal/mole,  jW  -  43.5  cal/ 

mole  “K,  and  c  ..W  38.5  cal/mole  “K;  lor  Nll.Cll).  (Hel.  12),  W  -  117.5  g/mole, 

p,  ff,  2  4  4 

=  513®  K  (a  solid  phase  transition),  but  since  heat  of  transition  data  and  heat  capacity 
data  appear  to  be  absent,  it  w,ll  be  assumed  that  Ah.  w  0 ,  c  i  **  2  **  0. 5  cal/g  ®K; 

and  =  300®  K  in  all  cases.  Utilizing  this  data  and  Eqs.  (2. 4)  through  (2. 6)  and  (2. 11) 
through  (2. 13),  it  is  possible  to  solve  for  as  a  function  of  Tp  for  all  three 

salts,  provided  AH  values  are  known. 


Using  the  values  of  AH  given  in  Eqs.  (2.  7)  through  (2.  9),  it  is  found  that  T  -  T 

P  s 

varies  from  1“  K  at  T  =  629®  K  to  14®  K  at  T  =  807“  K  for  NH  .Cl,  from  9®  K  at 
P  P  4 

T  =  453®  K  to  80®  K  at  T  =  573®  K  for  NH.NO.,,  and  from  3®  K  at  T  =  753®  K  to  40®  K 
P  p  4  3’  p 

at  Tp  =  903®  K  for  NH^CIO^.  Thus,  the  assumption  Tp  «  T^  is  approximately  valid  for 
NH^Cl,  invalid  for  NH^NO^,  and  poor  for  NH^CIO^.  It  is  interesting  that  the  variation 
in  Tp  -  Tg  is  such  that  the  experimental  checks  employed  indicate  that  the  approxi¬ 
mation  Tp«  Tg  is  valid;  e.g. ,  Tp  -  T^  is  relatively  small  at  the  melting  point  of 
NH  .NO„,  so  that  the  break  in  the  pyrolysis  curve  would  occur  at  a  value  of  T  near 

4  O  p 

the  melting  point. 


If  the  larger  AH  values  in  Eqs.  (2. 1)  through  (2.3)  are  used  instead  of  those  in 

Eqs.  (2.  7)  through  (2.  9),  then  the  calculated  values  of  T  -  T  are  roughly  doubled; 

P  ® 

it  is  clear  that  they  must  increase  since  the  term  mAH/W  enters  on  the  right-hand 

side  of  Eq.  (2. 11).  The  estimate  of  pD  Cp  used  in  this  calculation  should  not  be 

wrong  by  more  than  a  factor  of  5  and  therefore  cannot  substantially  modify  our 

conclusion.  Varying  i  does  not  alter  the  conclusion;  if  i  is  increased,  the  estimate 

of  T  -  T  increases  TEq.  (2. 10)1  ,  while  when  I  is  decreased,  free  molecule  heat 
p  s 

conchiction  sets  in  and  the  estimate  of  T  -  T  changes  very  little.  The  result  of 

P  s 

T  -  T  increasing  substantially  as  T  increases  for  NH.NO„  and  NH.CIO.  appears  to 

P  8  p  4  o  4  4 

be  inescapable.  A  recent  theoretical  analysis  by  Cantrell  (Ref.  13),  in  which  a  thin- 

film  lubrication  theory  with  heat  transfer  was  used  to  account  for  the  non-one- 

dimensional  flow  effects  in  pyrolysis  experiments,  also  supports  the  contention  that 

T  -T  is  not  small. 

P  8 
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The  main  effect  of  the  difference  T  -  T  is  to  modify  the  computed  activation  energy 

P  s 

for  the  interpretation  corresponding  to  Eqs.  (2.1)  through  (2.3)  and  to  modify  the 
AH  values  for  the  interpretation  corresponding  to  Eqs.  (2. 7)  through  (2.  9).  For 
NH^Cl,  NH^NOg,  and  NH^CIO^,  respectively,  the  new  activation  energies  are  roughly 
14  kcal/mole,  30  kcal/mole,  and  30  kcal/mole,  and  the  new  heats  of  sublimation  are 
very  roughly  14  kcal/mole,  25  kcal/mole,  25  kcal/mole. 

4.  FURTHER  COMMENTS  ON  AMMONIUM  PERCHLORATE 

Use  may  be  made  of  NH^CIO^  deflagration  e>?periment8  along  with  the  pyrolysis  experi¬ 
ments  to  provide  an  estimate  of  an  upper  bound  for  the  accomodation  coefficient, 
provided  the  surface  process  is  a  simple,  chainlike,  pure  sublimation  process 
[Eq.  (2. 9)]  in  both  pyrolysis  and  deflagration.  A  result  of  the  deflagration  analysis 
of  Johnson  and  Nachbar  (Ref.  14)  which  is  valid  for  this  purpose  is 

Yg  =  (Tf-Tg)/T*  (2.14) 

Certainly  <  Y  for  the  deflagrating  solid,  and  experimentally  T,  -  T„  « 116“  K, 

S  S 1 6  IS 

(T*  =  1415“  K).  This  yields  Y  ^  >  0.  082  at  1100“  K  and  100  atm.  Since^  «  1  in 

S>  6  P 

the  pyrolysis  experiments,  Eq.  (1.20)  implies 

L.e  ' 

for  simple  chainlike  surface  processes.  If  we  assume  l/ip  =  0  in  the  pyrolysis 

experiments,  then  Eqs.  (1. 16)  through  (1. 19)  imply  Y  =1  when  |  «  1,  whence 

s  p  p 

we  see 


s,e 


(2. 16) 


If  we  utilize  Eq.  (1.26),  with  AH  =  20  kcal/mole,  to  reduce  the  inequality  Y  >  0.082 

8 1  G 

at  1100“  K  and  100  atm  to  another  inequality  on  Y  at  900“  K  and  1  atm  (the  condi- 

S I  G 

tions  of  the  pyrolysis  experiment),  then  Eq.(2.16)  provides  an  inequality  for  [](4  /P)  +  ^  ] 
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Since  m  is  measured  in  the  pyrolysis  experiments,  and  /and  pD  are  easily  estimated 
(Section  3),  the  only  unknown  in  the  inequality  is  the  accommodation  coefficient  a. 

The  result  is  roughly  o  10  which  implies  that  the  pyrolysis  experiments  are  not 
near  surface  equilibrium.  All  of  these  conclusions  rely  strongly  upon  the  pyrolysis 
model  in  Part  One,  the  approximation  of  a  simple  chain-like  surface  sublimation 
process,  and  the  experimental  surface  temperature  estimates  in  deflagration. 
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